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Using the Bethe ansatz method, the zero frequency contribution (Drude weight) to the spin 
current correlations is analyzed for the easy plane antiferromagnetic Heisenberg model. The Drude 
weight is a monotonically decreasing function of temperature for all < A < 1, it approaches the 
zero temperature value with a power law and it appears to vanish for all finite temperatures at the 
isotropic A = 1 point. 
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The low frequency dynamics in one dimensional spin chains is a long standing problem. It has recently attracted 
a renewed interest, partly due to the fabrication of excellent quasi-one dimensional, spin 1/2 magnetic materials as 
SrzCuOz and CuGeO^. Detailed NMR experiments |l]] revealed an unusually high value of the spin diffusion constant 
and nearly ballistic behavior. 

The first issue on the question of spin diffusion is the zero frequency contribution (or Drude weight) to the dynamic 
spin current correlations at finite temperatures. If the Drude weight turns out to be finite, then the current correlations 
do not decay to zero at long times, implying ideal conducting behavior. If they decay to zero, the question still remains 
open whether they decay fast enough so that transport coefficients can be defined. Several numerical studies have 
been devoted to the analysis of the diffusive behavior in the Heisenberg model with suggestive but not conclusive 
results. 

In relation to this problem, it has been proposed that the integrability of the spin 1/2 Heisenberg model implies 
pathological spin dynamics and presumably the absence of spin diffusion m . A straightforward demonstration on the 
way in which conservation laws, characterizing integrable systems, might affect the long time dynamics was pointed 
out in reference H . There it was shown that in several quantum integrable models the uniform current correlations do 
not decay to zero at long times due to the overlap of the currents to conserved quantities. Unfortunately, this simple 
idea turned out to be insufficient for deciding about the decay of spin currents in the spin 1/2 Heisenberg model at 
zero magnetic field. 

On the other hand, a new method was proposed recently by Fujimoto and Kawakami |j] that allows the direct 
analytical evaluation of the Drude weight at finite temperatures. This procedure is based on the calculation of 
finite size corrections of the energy eigenvalues obtained by the Bethe ansatz method |J . The analysis starts from a 
convenient expression for the finite temperature Drude weight as the thermal average of curvatures of energy levels in a 
Hamiltonian subject to a fictitious flux coupled to the hopping or spin flipping term 1 1 Tjj| . Note that the anisotropic 



Heisenberg model is equivalent to the model of spinlcss fcrmions interacting with nearest neighbor interaction using 
the Jordan- Wigner transformation p^| . The direct analogy between charge and spin current correlations, suggests 
also the use of the name "Drude weight" in the context of spin correlations. 

In this work, we calculate the Drude weight for the antiferromagnetic Heisenberg model using the procedure proposed 
in references |]||. The formulation and notation of the thermodynamic Bethe ansatz equations by Takahashi and 
Suzuki will be closely followed. This construction is based on the string assumption for the excitations and it is 
particularly complex for arbitrary values of the anisotropy parameter A. The allowed type of strings are constrained 
by the nor maliz ability of the wavefunctions (TJ] . Therefore, for simplicity and without loss of generality, the analysis 
will be limited at A = cos(jr /v), ^=intcgcr, where only a finite number of string excitations is allowed. 

The results presented here are in good agreement with numerical results obtained by exact diagonalization of 
the Hamiltonian matrix on finite size lattices [0. They lend support both to the string construction and the novel 
procedure for calculating the Drude weight from finite size corrections to the Bethe ansatz eigenvalues. 

The XXZ anisotropic Heisenberg Hamiltonian for a chain of N sites with periodic boundary conditions S^ +1 = Sf 
is given by: 

N 

H = JJ2(S?S?+i + SfSf +1 + AS*S* +1 ) (1) 

i=l 

where Sf = \af , trf are the Pauli spin operators with components a — x,y,z at site i. The region < A < 1 is 
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parametrized by A = cos 0, 9 = ir/v, v = integer. The pseudomomenta k a and phase shifts <p a p characterizing the 
Bethe ansatz wavefunctions are expressed in terms of the rapidities x a : 

cot (y) = cot (^)tanh(^), 
cot(^) = cot(^)tanh(^M). (2) 
For M down spins and N-M up spins the energy E and momentum K are given by: 

M M 

E = J^(cosfc Q -A), K = Y,k a - (3) 

a— 1 a— 1 

Coupling the spin flipping term to a fictitious flux (p, the Hamiltonian becomes: 

N 1 

H = 3 H^°t°^ + h-c) + AS*Sf +1 ). (4) 

i=l 

The finite temperature Drude weight D can then be calculated by jll] : 
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1 ^ 1 &E n (cf>) i 
D =j\r2>5^H^°- (5) 



where E n are the eigenvalues of the Hamiltonian and p n the corresponding Boltzmann weights. Imposing periodic 
boundary conditions on the Bethe ansatz wavefunctions the following relations are obtained: 

sinh \ 9{x a +i) N _ e#JV ^ sinh \0{x a -x f3 + 2i) 12 M (6) 

X smh±6{x a -i) i l} 1 X sinh±d(x a -xp-2i) i ' ' '"' ' {) 

In the thermodynamic limit, the solutions of equations (||) are grouped into strings of order rij,j ~ l,...,u and 
parity vj = + or — . For 8 = -k /v the allowed strings are of order rij = j,j = 1, v — 1 and parity vj = + of the form: 

a£;+ = < + (« + 1 - M)i + 0{e~ 5N ); fc = l,2,...n, (7) 

and strings of order n v — 1 and parity = — of the form, 

x a - =x a +iu + 0(e- SN ), 5>0. (8) 

Multiplying the terms in equation ([|) corresponding to different members of a string and taking the logarithm we 
obtain: 

oo A/fc 

Nt j (xi)=2irli + ^2^2Q jk (xi-x>l)+n j ct>N-, a=l,2,...M j , (9) 

k=l /3=1 

I 3 a are integers (or half-integers) and M k is the number of strings of type fc, 

tj(x) = f&n^Vj), 
&jk(x) = f(x; \rij - n k \,VjV k ) + f(x;rij + n kl VjV k ) + 

Min(rij ,rifc ) — 1 

2 22 f{x;\nj -n k \ + 2i,VjV k ), 
i=i 

f(x;n,v) — 2utan~ 1 [(cot(n7r/2i/)) 1 ' tanh(7ra;/2j/)]. 

Following reference ||] the finite size corrections to the energy eigenvalues for a system of size N are calculated by 
introducing the function gij , g2j '■ 



x j N =xl + f + §t. (10! 
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where x : ' N (x J 00 ) are the rapidities for a system of size N(oo). Next, we expand equations (||) to orders of 1/N and 
in the thermodynamic limit introduce the densities of excitations pj and hole densities p^. The sums over the 
pseudomomenta are replaced by integrals over excitation densities plus boundary terms using the Euler-Maclaurin 
formula. 

To 0(1) we recover the integral equations for the excitation densities in the thermodynamic limit Ej ]: 

"j = \. (Pj + Pj) + X T i k * ■ ( U ) 

k 

* denotes the convolution a * b(x) — f^°° a(x — y)b(y)dy, Tjk(x) = (1 /2n)dO jk(%) / dx and djix) — (1 /2n)dtj (x) / dx . 
The sum over k is constrained to the allowed strings, given in our case by the equations and Xj = 1, j = 1, v— 1, 
K = -i 

To 0{l/N): 

*&j{Pj + Pj) = - E T ^ * Wk) + (12) 



To 0(1/N 2 ) 

^]92]{P] + Pj) + X T 3 k * (-92fePfc) : 
k 

l^9iM + p'D + T, T ^*^ikPk)} 



2dx 

k 



boundary terms (13) 



Minimizing the free energy we obtain the standard Bethe ansatz equations for the equilibrium densities rjj = p'j /pj 
at temperature T(f3 = 1/kbT): 

Inrjj = -2^sin(7r/z/)Ja J( 5 + ^A fc T^ * ln(l + (14) 

k 

These relations define the temperature dependent effective dispersions ej = (1/(3) ln(p^/Pj)- I n the string representa- 
tion the energy is given by: 

E = NY^ / dx{-2vs\xv{-)J aj {x))p 3 {x) (15) 

Expanding this expression for the energy we find that the first order correction in 1/N vanishes. So, the second 
derivative with respect to cj> of the second order correction gives us the final expression for the Drude weight: 

J ■> 

This expression is formally similar to the one obtained in reference Q for the Drude weight in the Hubbard model. 
It has an elegant interpretation by comparing it to the analogous expression for independent fermions. Taking the 
second derivative of the free energy with respect to the flux we find: 

% - E < "» > §f - < »» > c- < * »<fj? > 2 <"> 

where < n M > is the Fermi-Dirac distribution for particles with dispersion e M . Considering that the left hand side 
(the persistent current susceptibility ) vanishes in the thermodynamic limit for any finite temperature and that the 
first term in the right hand side is equal to 2ND, we find that: 
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D - h £ < > (I" < «M >)( af) 2 !^0 ( 18 ) 

Rewriting equation ( fL6| ) we arrive at a similar expression: 

^ = f_y x {p J+P )) < > (1- < »(^^) 2 (19) 

with < nj >= 1/(1 + e@ e i). So the Bethe ansatz expression for the Drude weight resembles that of independent 
fcrmion-like excitations. 

To obtain the distributions pj , Pj and , the coupled integral equations 
iteration. 

In Fig. 1, D is shown as a function of A for 2<^<16(0<A< 0.98) and different characteristic temperatures. 
The main result is that the Drude weight D is a monotonically decreasing function of A and temperature. At T = 0, 
D = f T(„-sa P"5[ - Most interestingly, D seems to vanish at all temperatures for A = 1. This result excludes an ideal 
conducting behavior for the isotropic Heisenberg model. Still, an anomalously slow long time decay of the current 
correlation functions could lead to pathological low frequency dynamics and non-diffusive behavior. Furthermore, the 
vanishing of the Drude weight at the isotropic point suggests that it remains zero at all temperatures in the region 
A > 1, the easy axis case (or insulating state in the fermionic model). This conclusion is in accord with the numerical 
results of reference [||. We should note that the numerical investigation close to the isotropic point is somewhat 
difficult as the number of equations to solve diverges. 

In the high temperature limit [fi — ► 0), D is proportional to (3. The constant of proportionality Cjj, equal to 
the long time asymptotic value of the current correlations [Q, is compared with results obtained in reference Q by 
exact diagonalization of the Hamiltonian on finite size lattices extrapolated to the infinite size limit. The quantitative 
agreement obtained lends support to the assumptions involved in the whole Bethe ansatz procedure for calculating 
thermodynamic properties and finite size corrections. 

The next observation is that the Drude weight approaches the zero temperature value with a power law of the form: 

D(T) = D(T = 0) - const.T", a = (20) 

To indicate this point, in Fig. 2, D(T — 0) — D(T) is shown for v — 3, ...,6 in a logarithmic plot along with lines of 
slope a. Note that the exponent a is half that for the low temperature spin susceptibility as obtained by Abelian 
bosonization [^6|. It is also consistent with the value a = 2 for free fermions {y = 2). 

The results presented above, concern only the zero frequency contribution to the spin current correlations. A 
reliable method for studying the low frequency behavior in intcgrable quantum many body systems (and the influence 
of non-integrable perturbations) remains a challenging problem. 



14),(0),(n2h are numerically solved by 
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FIG. 1. D(A) evaluated at the points v — 3, 16 and various temperatures. The continuous line is the high temperature 
proportionality constant Cjj = D/(3. The o indicate exact diagonalization results from reference (4). 

FIG. 2. D(T = 0) — D(T) at different temperatures (o's) in a logarithmic scale. The lines indicate slopes a = 2/{v — 1). 
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